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Abstract 

Using a geometric averaging procedure applied to a non-affine lin- 
ear connection, we prove that for a narrow one particle distribution 
function and in the ultra-relativistic limit, a bunch of charged point 
particles can be described by a Charged Cold Fluid Model, without 
additional hypothesis on the momentso 



1 Introduction 

Modeling the dynamics of a non-neutral plasma is an important problem 
in Plasma Physics and its applications in Beam Dynamics. Each bunch 
of a beam in an accelerator contains about 10 10 particles in a small region. 
Therefore, one looks for models such that: 1. Are simple enough to be useful 
in numerical simulations of the dynamics and 2. Contains the main features 
of the particle nature of the constitutive elements of the bunches. 
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The traditional approach has been to use fluid models. Since the bunches 
are formed by charged particles, a natural view is to interpret the fluid 
model as a macroscopic approximation of a kinetic model [1]. The resulting 
model is simple and still capable to retain some particle features coming 
from the underlying Kinetic Theory. However, these derivation is based on 
some assumptions that maybe are not accomplished in modern accelerators. 
These assumptions are usually on the high order moments of the distribution 
function. This seems to be a general feature of all the derivations of fluid 
models from Kinetic Theory [1]. 

We present in this note a new derivation of the cold fluid model from 
Kinetic Theory which uses only natural hypothesis, happening in current 
particle accelerator machines. The way to do this is the following. Firstly, 
we re-write the Lorentz force equation as an auto-parallel condition of a 
non-affine linear connection. Then, we use an averaging procedure described 
in [2] to average this connection. The resulting averaged connection is an 
afhne connection on the manifold M. After this, one compares the dynamics 
of both connections [3]. For narrow distributions which follow the Vlasov 
equation and in the ultra-relativistic dynamics, both dynamics are similar 
[3]. Therefore, we have two Kinetic Models (one based on the Vlasov equa- 
tion and the other on the Liouville equation associated with the averaged 
connection). It happens that under the same assumptions than for the par- 
ticle dynamics, the corresponding solutions of the Vlasov equation / and 
the averaged Vlasov equation / are similar [4]. One can also prove that the 
corresponding velocity vector fields are similar. Finally, the velocity field of 
the Liouville equation associated with the averaged dynamics is controlled 
by the diameter of the distribution and by the energy of the bunch [3] . Con- 
sidering all together, one has control of how good is the cold fluid model as 
an approximation of the Vlasov's model. 

Notation. The space-time manifold M is n-dimensional. The metric rj 
is the Minkowski metric, the potential A is a smooth 1-form. 

2 The Lorentz Connection 

The Lorentz force can be written in a covariant form in a general coordinate 
system in the following way: 
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i,j,k = 0,1, ...,n, 
(2.1) 
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where a : I — ► M is a solution curve for t £ I, are the coefficients of 

the Levi-Civita connection W of r/ and dA is the exterior derivative of the 1- 
form A (electromagnetic potential); here the parameter t is the proper-time 
of T] of the curve a. 

This system of second order differential equations are an auto-parallel 
condition of a non-affine linear Koszul connection [5] : 

L V & d = 0. (2.2) 

The connection coefficients are 

L r jk (x,y) = TV + - i J Fj^yV + F* 

t/ m i 
+ F 4 m (x) — (n,- fc — ; — ^ismy s y l ), (2-3) 

with Fjj := 9jylj — djAj and F* j = rj^F^j. This defines a non-affine con- 
nection on the pull-back bundle defined by the following diagram: 

f*TM^>TM (2.4) 



where 



N 2L_M, 

N:= [J {yG T X M, | ^(ar)yV' > 0}. 
xeM 



3 The averaged Lorentz dynamics 

The Lorentz connection L V is not an affine connection on M. One way to 
obtain an affine connection is to integrate over the support of the distri- 
bution / the connection coefficients [3]. This construction has a geometric 
interpretation [2], but is not relevant here. Applying this procedure to eq. 
(1.2), one obtains that the averaged connection < L V > has the following 
coefficients: 

< L r jk >= TV+r,- < - j^=, y m > v mk +F l k < - r^ y m > v mj )+ 

Wv{y,y) ^y/'niv^y) 

1 v m 1 

+F *-2(< m^^> )■ < ai > 
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Each of the integrations is equal to the y-integration along the fiber: 

< y i >■= ?/ L v / y i f(x,y)dfi, vol(E x )= f f{x,y)dn 
vol{2s x ) Js x J^ x 

and similarly for higher moments. Here I] is the unit hyperboloid, defined 
by 

£:= □ {y€ T x M\ Vij (x)yY = 1}. 
xeM 

T, x is the restriction of £ to T^M. 

Proposition 3.1 Let M and L V be as before. Assume that the support of 
the distribution f : XI — ► R is compact and denote by < L V > the averaged 
Lorentz connection. Then: 

1. The connection < L V > is an afftne, symmetric connection on M. 
Therefore, for any point x G M, there is a normal coordinate system 
such that the averaged coefficients are zero. 

2. To write down the form of < L V > we only need the first, second and 
third moments of the distribution function f(x,y). 

4 Comparison between L V and < L V > 

If on M there is a time-like vector field U normalized such that r/(U,U) = 1, 
one can define the Riemannian metric fj: 

fj(X, Y) := -n(X, Y) + 2 V (X, U) V (Y, U). (4.1) 

Using r) there is a scalar product on the vector space T^M defined by 
fj(x)ijdy l <S> dy 3 . The diameter of the distribution f x is 

a x := sup{dfj(yi,y2) 1 2/1,2/2 G sup(f x )}, a := sup{a x , x G M}. 

The Riemannian metric r) induces a distance function d^ on the manifold 
T X M. Then, we define a := sup{a x , x G M}. We choose as vector field U 
in the definition of the Riemannian metric (6.1) the following: 

U{X) = ^(x)</> y (x)<^>(x)> ' if KM < > (-) < V 3 > (-) > 

(4.2) 
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and if rjij(x) < y l > (x) < y° > (x) < 0. Given a continuous operator 
A x : T X M — ► T X M, its operator norm is 

\\AUx) := sup{ y£T x M\ {0} }. 

\\y\\v 

Let us denote by j(t) the gamma factor of the Lorentz transformation from 
the local frame defined by the vector field U to the laboratory frame, at 
some instance defined by the laboratory local time coordinate t. Denote 
by 9 2 (x) = ifix)- < y > 2 (x) and 9 2 (x) =< y > 2 (x) - y^ix). Here 
y{t) is the velocity tangent vector field along a solution of the Lorentz force 
equation and y{t) is spatial component of the tangent vector field along a 
solution of the averaged equation, with both solutions having the same initial 
conditions. The maximal values of this quantities on the compact space-time 
manifold are denoted by 6 2 and 6 2 . Then we can stay the following theorems 
[3]: 

Theorem 4.1 Let M be a semi-Randers space such that r\ is the Minkowski 
metric. Let us assume that: 

1. The auto-parallel curves of unit velocity of the connections L V and 
< L V > are defined for values of laboratory frame coordinate time at 
least t. 

2. The dynamics occurs in the ultra-relativistic limit, E(x) >> 1 for all 
x G M. 

3. The distribution function is narrow, oo > E(x) >> a for all x G M. 
4- It holds the following inequality: 

\e 2 -e 2 \ < i 

5. The support of the distribution function f is invariant under the flow 
of the Lorentz force equation 

Then, for the same arbitrary initial condition, the solutions of the equations 

L V x x = 0, < L V >, x = 

are such that: 

\\x(t)- x(t)\\ < (C(x)||F||(x) +C 2 (x)(l + B 2 {x)a))a 2 E- 2 {x)t 2 , (4.3) 

where C(x), C2(x) and ^(x) are functions on M and bounded by constants 
of order 1 . 
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Theorem 4.2 Under the same hypothesis as in theorem 4-1, the difference 
between the tangent vectors is given by 

xWW < {K(x)\\F\\(x) +K 2 2 (l + D 2 {x)a))a 2 E' 2 t. (4.4) 

with Ki and D 2 (x) functions bounded by constants of order 1 . 



5 Results on Kinetic Theory and Fluid Models 

In this section we overview our results from [4] on the approximation of the 
Vlasov Kinetic Model by Cold Fluid Model. 

In Kinetic Theory, the 1-particle distribution function follows the Liou- 
ville equation: 

xf(x,y) = o, (5.1) 

where \ is the Liouville vector field. In the case that the non-linear connec- 
tion is the Lorentz connection, we consider the Vlasov equation: 

L x f(x,y) = 0, (5.2) 

with the vector field 

L X = y k d k - L r k ij y i y j ^, i,j,k = 0,l,...,n. 

Similarly, for the averaged connection, the corresponding Liouville's equation 
is given by 

< L X > f(x,y) = 0, (5.3) 

with the vector field 

d_ 

Qyk ' 



< L X>-=y k d k -< L T k ij > 2/V70:, i,j,k = 0, 1, ...,n. 



This is a partial integro-differential equation where the averaged vector 
field < L x > is obtained using the distribution function associated with 
the Liouville vector field of the averaged connection < L V >. However, it 
is easier (because the existence of normal coordinates [5]) to extract some 
consequence from this equation. In particular, 

Proposition 5.1 Let f and f be solutions of the Liouville equations L x{f) = 
and < L x > (/) = such that a, the diameter of f is small; L \ and 

< L x > Q-fe the spray vector fields obtained from the connections L V and 

< L V >. Assume that a « E =< y° >j in the laboratory frame and the 
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averaged is performed using f . Then for the solutions of the Vlasov and the 
averaged Vlasov equation, 

\f(x(t),y(t)) - f(x(t),y(t))\ < C M ■ \x(t) - x(t)\ 

for some constant Cm depending on the manifold M. 

One also proves the following [4]: 

Proposition 5.2 With the above notation, the following relation holds: 

< L X >f = => < L V> V V= 0(a 2 ), (5.4) 
with a the diameter of the distribution f. 

The natural implication of the above results is the following [4] : 

Theorem 5.3 Let F be an external electromagnetic field and L V the asso- 
ciated Lorentz connection. Then the solutions of the Lorentz force equation 

L V ± x = (5.5) 

can be approximated by the integral curves of the normalized mean velocity 
field u of the distribution function f(x,y), solution of the associated Vlasov 
equation L \f = and the difference is controlled by polynomials on a at 
least of order 2. 

This theorem is the main result presented in this paper. It affirms the 
possibility of using the cold fluid model for narrow distributions and in the 
limit ultra-relativistic limit. 
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